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We study theoretically dynamical phases of vortices in superconducting films with arrays of ob- 
stacles. By performing a series of molecular dynamics simulations and analytical calculations, we 
demonstrate the existence of a phase of soliton-like vortex- density waves existing in a wide range 
of parameters. These waves are formed by a self- assembled phase separation process induced by 
strongly nonlinear density fluctuations of the moving vortex matter above a certain critical driving 
current. At high vortex concentrations, the waves move at an approximately current-independent 
speed resulting in a wide plateau in the volt age- current characteristics. At stronger drives, the 
vortex system enters into a fully jammed (zero- volt age) phase. By combining ac and dc drives, the 
interplay between the vortex- density- wave and jammed phases leads to the observation of negative 
absolute mobility of vortices, which induces the superconducting film into a negative resistance 
state. 



I. INTRODUCTION 

Driven many-particle systems are often affected by 
strongly nonlinear fluctuations which can lead to dra- 
matic phenomena such as self-organized density wavesi^, 
jamming-, and spontaneous segregation^. In inhomo- 
geneous type-II superconductors, nonlinear phenomena 
resulting from the motion of magnetic flux quanta (vor- 
tices) have been investigated, for instance, in the context 
of dendritic flux instabilities stemming from the break- 
down of the Bean critical stated—. 

In a superconductor, vortices are subjected to the 
Lorentz force induced by an applied current density J 
and to the pinning force produced by sample inhomo- 
geneities. A dissipative dynamical state of moving vor- 
tices is established at the critical current J c when the 
Lorentz force overcomes pinning £ Just above J c , mo- 
tion is usually plastic with vortices moving at differ- 
ent velocities, which results in a nonlinear velocity- force 
(voltage-current) characteristic. At a stronger drive, the 
pinning potential is washed out and the corresponding 
dynamical state, coined flux flow, is essentially linear. 
These dynamical regimes have been extensively inves- 
tigated in disordered superconductors^ and, more re- 
cently, in films with periodic arrays of nanoengineered 
pinning centers^— However, much less attention has 
been given to the dynamics of vortices in arrays of ob- 
stacles (or antipinning centers)^ Because in this config- 
uration vortices are not trapped individually, it renders 
a very small critical current, which, in principle, has lim- 
ited interest for applications. On the other hand a small 
J c allows for investigation of vortex motion in a much 
wider current range. Moreover, vortices may be forced 
into meandering paths inducing strong lateral fluctua- 
tions that can give rise to novel dynamical phases. 

In this paper, we demonstrate the breakdown of the 
flux-flow state in a superconducting film with an ar- 
ray of obstacles into a vortex- density wave state and, 
subsequently, into a fully jammed phase. Evidence of 
density waves in vortex systems have previously been 



found in the context of current-depaired vortices and an- 
tivortices in clean superconductors 14 and turbulence in 
superfluids^. In contrast, the phenomena we present 
here are collective states of the moving vortex matter 
resulting from the interplay of vortex- vortex interactions 
and the friction induced by the obstacles. As we shall see 
in this paper, these competing interactions lead to a crit- 
ical profile of the density waves where distinct dynamical 
phases, as well as a jammed phase, coexist. 

The paper is organized as follows. In Sec [III we gi ye 
the details of our model and numerical procedure. In 
Sec. IIII( we present the main results of our molecular 
dynamics simulations and discuss the dynamical phases. 
An analytical model is proposed in Sec. [TV] to explain 
the conditions for formation and stabilization of vortex- 
density waves and how this leads to a constant voltage 
regime. Sec. IIVI is devoted to illustrating an application 
of the interplay between the moving and jammed phases, 
namely, the negative absolute resistance effect. Finally, 
our main findings and final remarks are summarized in 
Sec. EH 

II. MODEL AND NUMERICAL DETAILS 

We consider a set of N vortices generated in the film by 
a perpendicular magnetic field B = Bz. The dynamics of 
a vortex i is modeled by the Bardeen-Stephen equation, 

VVi = F-ViU b -J2ViU ij: (1) 

where r] is the viscous drag coefficient, F is the Lorentz 
force induced by the applied current, and Uij is the vor- 
tex pair potential, modeled here as Uij = eKofaj/A) 
(e = 0o / '(47r/ioA) and </>o is the flux quantum). A is 
the effective penetration depth measuring the range of 
vortex- vortex interactions. For A much larger than all 
length scales of the system Uij asymptotically reduces 
to a logarithmic potential. In this investigation, we an- 
alyze the effect of short and long range interactions by 
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FIG. 1. (Color online) Foreground: unit cell of the triangular 
array of ferromagnetic barriers. Background: contour plot 
of the potential Ub(x,y) generated by the array. Each cross 
comprises two legs perpendicular to one another, of length 
0.5a, width 0.04a, and height 0.025a. To prevent proximity 
effects, the ferromagnets are assumed to be separated from 
the superconductor by a 0.0125a thick insulating layer. 



conducting simulations for 0.5 < A < oo. Hereafter we 
adopt the following units: a (obstacle lattice constant) 
for length, to = e/(r]a 2 ) for time, and e for energy. 

Ub{x, y) is the potential resulting from a triangular ar- 
ray of cross-shaped obstacles (Fig. [I]). To be specific we 
have chosen as such obstacles small ferromagnets (FMs), 
each with a uniform permanent magnetization M an- 
tiparallel to B. This choice is motivated by the well- 
known properties of FMs to either attract or repel vor- 
tices, depending on their magnetic orientation^. Other 
possible choices include higher- T c superconducting dots 17 
and pillars^. Each FM repels a vortex j via a potential 
given by U kj = - f M • bjd 3 r k , where f k (x k ,y k ) is 
the position of the volume element of the micromagnet k 
and bj is the local flux density of vortex j at f k calculated 
within the London approximatio n 18 ! 19 . In the results pre- 
sented here we used M = 320(/>o/Mo& 2 (= 0. 8^2/^0, 
assuming a coherence length £ = 0.02). For this value 
the ferromagnets do not induce any vortex- ant ivortex 
pairs 2 ^. Ub(x, y) results from the superposition of all FMs 
in the array (Fig. [I]). It introduces a current-induced 
trapping mechanism between two critical drive values: 
F c i = 3.5, at which a vortex can get trapped at a cross 
corner, and F C 2 = 12.0, above which a vortex surmounts 
the barrier. These properties provide the main ingre- 
dient for the phenomena we investigate here, namely, a 
non-linear, drive- dependent friction. 

Eqs. [1] are solved numerically via molecular dynam- 
ics (MD) simulations in a cell of size L x x L y with pe- 
riodic boundary conditions. The values of L x and L y 
were chosen after a careful finite size analysis. We veri- 
fied that the transitions between the observed dynamical 
phases become unaltered for system sizes typically larger 
or equal to 24x48 or 12x96 FMs. Therefore, N typically 
ranges up to a few thousands of vortices. The simula- 
tion procedure is as follows. First, the vortex system is 
equilibrated at zero current via a standard simulated an- 
nealing schem o 2Q i 21 . Then, the driving force F yi applied 
along the y axis is slowly increased. For each F yi time 
series and averages are calculated on an interval At of 



typically 10 6 time steps after a stationary state has set- 
tled. The main physical quantities are the center-of-mass 
velocity v y = ^ . v { • y of the vortex array along the 
drive direction and its time average. 

III. DYNAMICAL PHASES 
A. General picture 

We conducted an extensive series of simulations of the 
model described above for densities ranging from n = 0.1 
to 2.0 vortices per obstacle and different values of A. 
Fig. [2]- (a) presents (v y )-F y curves for a few n values. 
The different dynamical behaviors are arranged in the 
diagram of Fig. [2]-(b). In general, for F < F&, all vor- 
tices meander through the potential channels leading to 
a linear (flux-flow) phase, which we call phase I. This 
phase has an effective viscous drag coefficient, defined 
as ?7eff = F/(v), approximately force independent and 
smaller than 77. This is consistent with the fact that the 
meandering motion of vortices dissipates more energy as 
compared to the case of conventional flux flow. 

Above F c i, we observed a remarkable plateau in 
(v y )(F y ) for long-ranged vortex- vortex interactions. This 
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FIG. 2. (Color online) (a) Vortex mean velocity (v y ) as a 
function of the driving force F y for logarithmic (A = 00) 
vortex-vortex interaction and occupation numbers ranging 
from n — 0.25 to 1.5. The roman numbers correspond to 
the phases depicted in the F y -n diagram (b) and explained 
in the text, (c) (v y ) versus F y for n = 0.25 and different A 
values: A = 0.5 to 4 (symbols) and A = 00 (line). 
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FIG. 3. (Color online) Snapshots of the flux density distribution b(x,y) in a 12x144 array of cross barriers for total vortex 
density n — 0.25, A = 4.0 [(a)-(e)] and 1.0 [(f)-(j)], and different force values spanning phase I [(a) and (f)], II [(b)-(d) and 
(g)-(i)], and III [(e) and (j)]. The drive direction (y axis) points left to right. The curves on the right depict the flux density 
b(y) (in units of 0o per unit cell) integrated over the sample width. The shaded areas represent the shape of the density wave 
as predicted by the mean-field model presented in Sec. IIV1 The plots on the left represent a zoom-in of the regions indicated 
by the squares showing the positions of vortices (dots) and barriers (crosses). 



constant velocity state, phase II, persists in a wide range 
of force values and ends drastically in a fully jammed 
((v y )=0) state for n < 0.75, phase III, or in a mov- 
ing linear phase IV for larger n. At even higher drives 
[F > F C 2), all vortices surmount the barriers leading to a 
new flux- flow regime (phase V). Phase II dominates the 
diagram in the range of parameters studied. To check the 
stability of this phase with respect to the kind of vortex- 
vortex interaction, we fixed n = 0.25 and run simulations 
for A values down to 0.5. As observed in Fig. [2]-(c), the 
plateau shrinks as A decreases and is substituted by a 
smooth decrease of (v y ) for A < 1.0. As we shall see be- 
low, this region of decreasing (v y ) have similar properties 
to the plateau region of the A > 2 systems. 



B. Vortex-density waves 

Here we analyze phase II in more detail. For simplicity 
we focus only on the vortex density n = 0.25. In Fig. [3j 
we plot snapshots of the flux density distribution b(x, y) 
for A = 4 (long range) and A = 1 (short range). The cen- 
tral panels correspond to density plots of y) whereas 
the panels on the right are flux profiles averaged over the 
x direction, i.e., b(y) = j dxb(x,y) / L x . The left panels 



depict the snapshots of vortex positions in the areas in- 
dicated. These plots reveal a spontaneous breakdown of 
the vortex distribution when going from phase I [panels 
(a) and (f)] to phase II [(b)-(d), for A = 4.0, and (g)-(i), 
for A = 1.0]. The resulting patterns comprise stripes of 
high flux density oriented perpendicularly to the drive 
direction and traveling at constant speed as soliton-like 
vortex- density waves (VDW). 

For long range interactions, the waves have a peculiar 
triangular flux profile, which changes shape as the force 
is increased, and are intercalated by narrow regions of 
sparsely distributed vortices trapped by the cross cor- 
ners. By carefully analyzing the dynamics in different 
regions of the sample, we found out that the left side of 
the wave is characterized by meandering motion of vor- 
tices, similar to the dynamics at F < F& (phase I), while 
on the right side vortices will rather assist each other to 
jump over the barrier, which is the same dynamics found 
at F > F C 2 (phase V). Therefore, a single density wave 
corresponds to a moving region where two distinct dy- 
namical phases coexist. Such phase separation was also 
observed for short range interactions, but not so clearly 
because in this case a density wave spans only a few lat- 
tice spacings of the FM array. In addition, the pinned 
vortex regions for small A are considerably larger and in- 
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FIG. 4. Power spectra of density fluctuations at different 
drives for n = 0.25 and A = 1.0 (a) and A = 4.0 (b). The 
curves are displaced from each other by two decades for bet- 
ter visualization. In both panels 1/f 2 curves (gray line) are 
plotted as guides to the eye (see text). Insets: driving force 
dependence of the mean vortex velocity and the density wave 
velocity for A = 1.0 and 4.0. 



creases with force, which is accompanied by a reduction 
of the mean vortex speed. A detailed analysis of the mor- 
phology of the VDWs and its dependence on the driving 
force and other parameters is given in Sec. HVl 

Although the vortex wave patterns observed in our 
simulations are not space periodic, they do present time 
order. A convenient way to search for time correla- 
tions in the vortex density n(t) is by studying the power 
spectrum of its fluctuations, S(f) = (\n(f)\ 2 ), with 
n(f) = Jdt n(t)e~* 27r ^ £ £ We perform this task by record- 
ing the density n(t) integrated over a 12 x 12 region of 
the film every other 10 time steps during a series of 10 6 
time steps. S(f) is then estimated by averaging \n(f)\ 2 
over N s = 10 segments of the time series. 

Plots of typical £(/), for both short (A = 1.0) and long 
(A = 4.0) range cases, are depicted in Fig. 2J As it is 
clear in this figure, S(f) has a long 1/f 2 tail for all drive 
intensities and interaction ranges studied. This high fre- 
quency behavior is associated with Brownian noise re- 
sulting mainly from the meandering dynamics, where a 
vortex colliding with a barrier has to choose randomly 
whether to dodge left or right. The low frequency re- 



gion of S(f), however, changes dramatically from a fiat 
behavior in phase I to a series of strong peaks in phase 
II. These peaks correspond to the harmonics associated 
with the coherent passage of high flux density domains at 
constant speed. Therefore, the power spectra provide an 
efficient means of determining the density-wave velocity: 
v w — f\L y , where f\ is the first harmonic frequency. 

We have calculated the wave speed for different force 
values both for A = 1.0 and 4.0. For better accuracy, we 
calculate v w from a (sharper) higher harmonic f s = 5/1 
(typically, s = 3 or 4). The results are plotted in the 
insets of Fig. 3J Remarkably, v w is always larger than 
the mean vortex speed for both kinds of interactions. 
For short-range interaction, v w increases with the driving 
force, whereas (v) decreases. Such high speeds of the 
density waves as compared to (v) can be explained by 
the presence of slow moving vortices between adjacent 
waves. In fact, because vortices in a density valley are 
diluted, they can easily be trapped at barrier corners 
(thus contributing to decrease (v)) until being catch by 
a wave front. In opposition, vortices in the density peaks 
assist each other to overcome the traps. Put in another 
way, the VDWs observed here are anti- jamming waves 
self-assembled in a way as to avoid a full stop of the 
vortex flow. 

Density waves have been observed in other classical 
many-particle systems, for instance, car tram o 22 i 23 , gran- 
ular flo w 24 i 25 and galaxies 2 ^, and their origin is still un- 
clear. In contrast to the anti-jamming waves observed 
in our simulations, the density peaks in these systems 
are jamming domains intercalated by free motion re- 
gions. Accordingly, their wave velocities are usually 
smaller than the mean flow velocity and can even be- 
come negative^. However, in all these systems, ours in- 
cluded, density waves can be understood as a dynam- 
ical, self-organized phase separation, similar to equilib- 
rium pattern formation in systems with competing in- 
teractions^. For instance, in the case of granular media 
flowing through narrow tubes 2 ^ density waves form spon- 
taneously and generate a spectrum of density fluctuations 
very similar to those presented here in Fig. 2J in spite of 
the very different types of interaction in these systems 
and different nature of the waves. 



IV. MEAN-FIELD APPROACH 

To understand how the vortex-density waves are 
formed and why they lead to a constant voltage regime 
for large A, we propose an analytical model based on a 
mean-field analysis of Eq. [TJ We start from the Fokker- 
Planck equation for the vortex distribution function at 
zero temperature and coarse-grain it over scales larger 
than A. This procedure leads to the following equation 
for the vortex density p'£^ 



-V- 



p(P- gVp + Fp) 



(2) 
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—gVp is the lowest order (long-wavelength) contribution 
of vortex- vortex interactions, with g = — ^ J cPr r-VUij. 
This indicates that, as expected, a local flux gradient 
introduces an additional term to the local driving force, 
that is Fd = F — gV p. Fp is the coarse-grained force 
produced by the pattern, which can be understood as a 
drive-dependent friction force. 

To account for a dynamical trapping mechanism sim- 
ilar to that induced by the cross-shaped barriers, we as- 
sume that Fp represents a static friction in the range 
Fd < Fd < F C 2, for which no motion occurs. For Fd < 
F c i and Fd > F c2l motion takes place and, accordingly, 
Fp must represent the corresponding kinetic frictions F^i 
and Fjz2, respectively. These properties can be arranged 
by modeling the friction force as Fp = —yFp(Fd), with: 

{F kl , HF d <F cl 
F P =<F k2 , HF d >F c2 (3) 
[F d , ]£F cl <F d <F c2 . 

We proceed by assuming a one-dimensional density 
wave is formed at a given F, with F& < F < F c2 , and 
search for the conditions for such a wave to exist with a 
fixed shape and a constant velocity v w directed parallel 
to F = yF. In this case, we can transform into a refer- 
ence frame moving with the wave, (y f = y — v w t, t' = £), 
and rewrite Eq. [2] as 

rjv w = F - gd y >p- F P: (4) 

where p is assumed to be a function of y' only and d y > p 
is the density gradient in the moving frame. To deter- 
mine the shape of the wave, it is convenient to choose 
the position of the wave maximum as the origin of y' 
and analyze separately the regions I = {y f < 0} and 
r = {y' > 0}. Further, for simplicity, we assume the 
wave is single-peaked and define G r ^ = \d y tp\ r (i), such 
that F d = F - gGi , for y' < 0, and F d = F L + gG r , for 

y'>o. 

With the help of Eq. O we can state the following 
conditions of motion of the wave: 

gGi>F-F cU (5) 
gG r > F c2 - F, (6) 

where we used the fact that, for any F G [F c i,F c2 ], the 
positive slope on the right side of the wave rules out the 
possibility Fd > F c2l whereas the negative slope on the 
right rules out Fd < F& . This naturally accounts for the 
dynamical phase separation observed in our MD simula- 
tions and discussed in Sec. IIIIB1 

A direct consequence of the above condition is that if 
the vortex distribution peaks up at a certain point, the 
corresponding bump can only turn into a traveling wave 
if the gradients in both sides of the bump fulfill Eqs. [5] 
and[6j Subcritical bumps, i.e. those with gG\ < F — F& 
and gG r < F c2 — F, would simply be trapped by the 
pattern. On the other hand, an increase in G r above its 



critical value leads to an increase in the right side velocity 
while an increase in G\ above F c2 — F reduces velocity on 
the left side. Therefore, if the vortex distribution peaks 
up too sharply the wave will be unstable and will rapidly 
decay to a profile with gentler slopes so that the velocities 
on the left and right sides of the wave come to a common 
value. It is clear then that the condition for the stability 
of a wave traveling at constant speed and fixed shape 
is that the density gradient must be everywhere critical, 
that is 

Gi = (F L -F cl )/g, (7) 
G r = (F c2 -F L )/g. (8) 

This result is analogous to the build up of the Bean crit- 
ical state in the magnetization of hard type-II supercon- 
ductors, where the flux gradient (induced by the external 
flux penetrating the material) exactly balances the crit- 
ical force induced by pinning centers. Here the critical 
forces are exactly balanced by the overall driving force, 
which accounts for the applied Lorentz force and the flux 
gradient. 

Eqs. and [8] describe a triangular shape for the density 
wave which switches its left (right) side slope from gentle 
(steep) to steep (gentle) as the Lorentz force is increased 
from F c i towards F c2 . To contrast this result with our 
simulation data, we take F& = 3.5 and F c2 = 12, ob- 
tained from the (v)(F) characteristic for A = 4 and 
n = 0.25, and calculate the parameter g for Bessel-like 
vortex- vortex interactions, which gives g = 2irh?e. Then, 
we calculate the critical gradients from Eqs. [71 and [8] and 
compare with the flux profiles obtained from the MD sim- 
ulations for F = 5.0, 7.0, and 9.0. The results, presented 
in Fig.[3](b)-(d), demonstrate a good agreement between 
the numerics and the mean-field model. 

To estimate the wave velocity, we first impose that the 
velocity in the left side of the wave (yi = v w ) must be con- 
sistent with the linear, meandering dynamics expected to 
occur for Fd < F c \. For this kind of motion, an increase 
in the drive leads to a proportional increase in the veloc- 
ity, i.e., v w = Fd/r] e ff. Then, we notice that, from Eq. [71 
Fd = F — gGi = F c \. Hence, the wave velocity does not 
depend on F: 

Vw = ^cl/^eff- (9) 

This is in excellent agreement with the MD results for 
large A if one assumes that all vortices participate on 
a wave, that is (v) = v w . In fact, as revealed by 
the simulations, there are narrow, low-density regions of 
trapped vortices coexisting with the dynamical phases 
for n = 0.25 and large A. For that reason, v w is typically 
slightly larger than (v). 

It is worth noticing that in the limit F& -^0a moving 
wave solution is not possible (at least in the framework 
of the present model). Therefore, a condition for the oc- 
currence of vortex-density waves is that vortices interact 
with a potential which is able to trap them only when the 
applied current surpass a certain critical value. This is 
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probably the reason why vortex-density waves have never 
been detected or predicted in, for instance, arrays of pin- 
ning centers. In these systems, in general, there is no 
moving phase preceding the pinned phase. 

It is also important to emphasize that the mean-field 
approximation (Eq. [2]) is strictly valid for vortex densi- 
ties n ^> A -2 . In the simulations performed for A = 4a 
and n = 0.25^ a~ 2 we have n = 3. 46 A -2 , which can 
be considered inside the range of validity of the model. 
Indeed, the main results obtained in the MD calcula- 
tions for these parameters were correctly reproduced by 
the mean-field approach. For A = a, however, the vor- 
tex distribution is too diluted. In addition, in order to 
treat the interaction with the barriers macroscopically, A 
should span several lattice spacings. For these reasons, 
the analytical results presented here are not strictly valid 
for short range interactions. Notwithstanding, our sim- 
ple analysis ignores the details of the barriers, taking into 
account only its main macroscopic properties. Therefore, 
it can be applied to other barrier configurations, as long 
as their coarse-grained properties can be expressed in a 
way similar to Eq. [3j 



V. NEGATIVE ABSOLUTE RESISTANCE 

The dynamical transition between the VDW and jam- 
ming phases opens the possibility for an interesting appli- 
cation: the construction of an active device (made of su- 
perconducting material) exhibiting negative absolute re- 
sistance (NAR). The working principle of this device is 
based on the negative absolute mobility (NAM) effect, 
which corresponds to motion in a direction that opposes 
the driving force, irrespective of the drive direction. This 
phenomenon was predicted to occur in the transport of 
a single classical Brownian particle through a symmet- 
ric, periodic substrate 2 ^ and was experimentally corrob- 
orated in a system of colloidal spheres in a microfluidic 
device^ and subsequently in a Josephson junction^. In 
the later, the phase dynamics, which can be mapped into 
the problem of a single Brownian particle, lead to nega- 
tive resistance. 

Quite generally, the main ingredients for the observa- 
tion of NAM are: (i) a medium that allows easy motion 
of the particles at low drives and strongly suppresses mo- 
bility at high drives; (ii) a fluctuating (nonequilibrium) 
force, for instance an ac excitation, superposed to the 
dc drive. The first requirement implies necessarily that, 
for a certain force range, the mean velocity of the driven 
system must decrease as the driving force increases, that 
is, negative differential mobility. This phenomenon have 
been experimentally demonstrated for vortices driven in 
superconducting films with periodic arrays of pinning 
centers at magnetic fields close to one flux quantum per 
pinning site (Bi)^ Previous numerical calculations have 
pointed out that such negative differential mobility of 
vortices near B\ is a result of a transition from a highly- 
dissipative, disordered regime to a filamentary state with 
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FIG. 5. Vortex mean velocity as a function of the dc force 
Fdc in the presence of a sinusoidal excitation with frequency 
uo — 0.01 and amplitudes A — 11.25,10.25, 7.5 and 5.5 corre- 
sponding to A = 4,2, 1,0.5, respectivelly. All curves exhibit 
absolute negative mobility of vortices at small Fdc- 



somewhat smaller but non- vanishing mean velocities J^^ 2 - 
As shown in Fig. [2j the VI characteristics of the vortex 
system under investigation here do present negative dif- 
ferential resistance, in a broad field range, as a result of 
the transition from the VDW to the jammed state. Here, 
however, the voltage drop corresponds to a strong sup- 
pression of vortex mobility, which, as discussed below, is 
a crucial feature for the observation of NAM. 

Here we demonstrate how the NAM effect can be ap- 
plied to Abrikosov vortices in patterned superconduc- 
tors in order to achieve negative resistance. We subject 
the vortices to a Lorentz force F y = F& c + F ac (t), with 
Fac(t) = A sm(ut) and an ac amplitude A fixed at a value 
chosen just above the VDW-jamming (II-III) transition, 
and measure (v y ), which is proportional to the dc volt- 
age, as a function of the dc drive F^ c , proportional to the 
dc current. The ac term of F y provide the extra energy 
source required for the observation of NAM. The results 
are shown in Fig. [5] for several A values. At zero dc 
current (F^ c = 0) the vortex dynamics cycle symmetri- 
cally through phases I, II, III, and back, producing zero 
dc voltage. However, for small positive (negative) F^ c , 
vortices are most of the time in the pinned phase III dur- 
ing the positive (negative) half-cycle of F ac and always in 
one of the moving phases I or II in the negative (positive) 
half-cycle. This leads to a net motion contrary to the dc 
force, that is, the mean electric field induced by vortex 
motion points antiparallel to the applied dc current, thus 
generating the negative absolute resistance effect. 



VI. CONCLUSIONS 

In conclusion, we studied novel dynamical phases 
of vortices in a patterned superconducting film. By 
molecular dynamics calculations we demonstrated, for 
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a wide range of parameters, the existence of vortex- 
density waves propagating at a constant speed v w larger 
than the mean vortex velocity (v) and essentially force- 
independent for long-range vortex- vortex interactions. 
The waves consist of well-defined regions of coexisting 
moving phases, where vortices assist each other to either 
dodge or overcome the barriers, and are intercalated by 
regions of trapped vortices. 

Our mean- field analysis revealed that these waves stem 
from fluctuations in the vortex distribution induced by 
a highly nonlinear friction force. A density bump pro- 
duced by such fluctuations turns into a stable moving 
wave when it reaches a certain critical profile. It is pre- 
cisely this critical shape of the waves that results in a 
force independent wave velocity, thereby naturally ac- 
counting for the voltage plateau observed in the MD sim- 
ulations. 

We have also demonstrated the feasibility of a super- 
conducting device which exhibits negative absolute resis- 
tance. This is achieved by exploring a combination of 
ac and dc excitations in such a way as to conveniently 
switch the system between the moving and fully-jammed 
phases. This prediction could be promptly tested us- 
ing conventional transport measurement techniques on 



a nanostructured sample with a pattern similar to that 
proposed here. 

Finally, it is worth pointing out that, given the gener- 
ality of our mean field model, the main results predicted 
here can also be applied to other systems of interacting 
particles, such as colloids and pedestrians. An advantage 
of the vortex system proposed here is that interactions 
can easily be tuned. In typical nanostructured supercon- 
ducting films, A can be varied from a fraction to sev- 
eral lattice spacings by controlling the film temperature 
near T ^ 21 i 33 , thereby allowing for experimentally access- 
ing our predictions in both long and short range cases. 
In additional modern imaging techniques could be used 
to identify the vortex density waves. 
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